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The Decay of Isotropic Turbulence 
in a Rapidly Rotating Frame 


By C. G. SPEZIALE 1 , N. N. MANSOUR 2 , and R. S. ROGALLO 2 


A direct numerical simulation of the decay of initially isotropic turbulence in a 
rapidly rotating frame was conducted. This 128x128x128 simulation was completed 
for a Reynolds number Re\ — 15.3 and a Rossby number Ro\ = 0.07 based on 
the initial turbulent kinetic energy and Taylor microscale. The numerical results 
indicate that the turbulence remains essentially isotropic during the major part of 
the decay (i.e., beyond the point where the turbulent kinetic energy has decayed 
to less than 10% of its initial value). The rapid rotation has the primary effect of 
shutting off the energy transfer so that the turbulence dissipation (and hence the 
rate of decay of the turbulent kinetic energy) is substantially reduced. Consequently, 
the anisotropy tensor remains essentially unchanged while the energy spectrum 
undergoes a nearly linear viscous decay — the same results that are predicted by 
Rapid Distortion Theory which is only formally valid for much shorter elapsed times. 
Surprisingly, no Taylor-Proudman reorganization of the flow to a two-dimensional 
state is observed. The implications that these results have on turbulence modeling 
are discussed briefly along with prospective future research. 


1. Results 

The research conducted this summer at the CTR concentrated on the develop- 
ment of improved Reynolds stress models for the description of rotating turbulent 
flows. It is envisioned that such models could also have important applications in 
the description of curved turbulent flows as a result of the analogy that can quite 
often be drawn between rotation and curvature. 

In order to gain insight into the effects of rotation, a direct numerical simula- 
tion of decaying isotropic turbulence in a rapidly rotating frame was conducted. A 
Reynolds number of Re\ = 15.3 based on Taylor microscale and a Rossby number 
of Ro\ = 0.07 were considered (this Rossby number is more than an order of magni- 
tude smaller than those which were considered previously). This direct simulation 
yielded some surprising results. As has been shown in previous numerical simu- 
lations and experiments (see Bardina, Ferziger and Rogallo 1985, and Wigeland 
and Nagib 1978), the turbulence remained isotropic after the rotation was imposed. 
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FIGURE 1. Transfer spectra as function of k i2 = ykj + and k 3 . a) Initial 
distribution (isotropic decay), b) Shortly after the rotation is started (ito / <?o )• 


The rotation killed the energy transfer (through the generation of inertial waves; see 
Fig. 1) and the turbulence underwent a pure viscous decay as would be predicted 
by Rapid Distortion Theory (RDT). More precisely, the energy spectrum, E(k,t), 
decayed in time in good agreement with the formula, 

E(k,t) = E(k,t 0 )ex p [-2 uk 2 (t - t 0 )] (1) 

where u is the kinematic viscosity (see Fig. 2). Equation 1 is obtained from RDT 
for this problem. The surprising finding was that Eq. (1) remained an excellent 
approximation even after the turbulent kinetic energy had decayed to only 10% of 
its initial value. As a result of the energy transfer being suppressed, the turbulence 
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Figure 2. Energy spectra, E uu , E vv and E ww (g 2 /?o = 0.03). 


decayed slower in the rotating frame (see Fig. 3). Although the integral length 
scales showed the development of mild anisotropies, there was no discernable Taylor- 
Proudman reorganization to a two-dimensional flow. The tensor, 


Aij 




(2) 


which is normalized by the dissipation rate e, remained isotropic (under a complete 
Taylor-Proudman reorganization, the velocity gradient along the axis of rotation 
u i3 i-» 0 as the rotation rate fl oo and, hence A 33 << A u , A-^). It can be shown 
that the RDT solution does not undergo a Taylor-Proudman reorganization since 
the Fourier transform of the velocity, 


U{(k,t) oc. A(k) exp(ia(k)Clt) 
and, hence, in the limit as fl i— > oo, 


fl dt 


0 ( 1 ) 


(3) 

(4) 


where u is the vorticity vector. 

In order for the Taylor-Proudman theorem to apply, (l/Sl)du>i/dt must vanish as 
fl oo. Since RDT becomes a better approximation for longer instants of time 
as fl gets larger (for a given turbulence level), it appears that no Taylor-Proudman 
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FIGURE 3. Effect of rotation on the decay of the turbulent kinetic energy. 


reorganization will occur in a rapidly rotating isotropic turbulence. Previous in- 
vestigators (c./., Bardina, Ferziger and Rogallo 1985) had speculated that such a 
reorganization to a two-dimensional flow could occur. 

Since the results of these direct simulations on rotating isotropic turbulence 
clearly demonstrate (in support of Wigeland and Nagib 1978, and Bardina, Ferziger 
and Rogallo 1985) that there is a reduction in the dissipation rate with increasing 
rotation rate, it is clear that modifications need to be made in the dissipation rate 
equation. All of the commonly used dissipation rate equations (c./., Launder, Reece 
and Rodi 1975) predict that for a given mean flow, a system rotation has no effect 
on the evolution of the dissipation rate in contradiction of experimental and numer- 
ical simulation data. A recently proposed model by Bardina, Ferziger and Rogallo 
(1985) given by, 

€ = —C\—z — C 2 fie (5) 

r 

for an isotropic turbulence (where e is the dissipation rate, q 2 is the trace of the 
Reynolds stress tensor, and C 2 = 11/3 and C\ — 0.15 are empirical constants) 
was tested. It was found that this model, which compared favorably with the data 
of Wigeland and Nagib (1978) performed very poorly at the rapid rotation rates 
considered herein (see Fig. 4). Consequently, it appears that the dependence of the 
dissipation rate equation on Cl is unlike in Eq. (5). Furthermore, it is not clear 
at this time how such a modified dissipation rate equation could be generalized to 
anisotropic turbulent flows. Any smooth function of the invariants, 
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Figure 4. Comparison of the decay of the turbulent kinetic energy with the 
model of Bardina, Ferziger & Rogallo (1985). 


where, 

toij = 2 (Ui,j ~ Uj t i ) , Sij = - ( Ui t j + Uj t i) (7) 

Ui is the mean velocity, and Tjj is the Reynolds stress tensor, reduce to functions 
of ft in the limit of rotating isotropic turbulence (a model containg tojjSij was 
proposed by Pope 1978). Work on the development of a more general dissipation 
rate equation which can account for rotation in inhomogeneous turbulent flows was 
begun during the summer program and will continue in the future. 

2. Future work 


The time development of the energy spectrum for an isotropic turbulence is given 

by, 

E(k,i) = T{k,t)-2vk 2 E(k,t) (8) 

where T(k,t) is the energy transfer. The equation for the dissipation rate, c, can 
be derived by multiplying Eq.(8) by 2vk 2 dk and integrating over all wave numbers, 


e = 2u ( k 2 T(k,t)dk - f k 4 E{k,t)dk 
Jo Jo 


( 9 ) 


The commonly used model for the right-hand-side of Eq. (9) is given as follows: 


2 vjf k 2 T(k t t)dk — 4i/ 2 f k 4 E(k,t)dk 



( 10 ) 
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where q 2 is twice the turbulent kinetic energy. Any modification to the above model 
that takes into account the effect of rotation, has to reflect the fact that the first 
term in Eq. (10) vanishes for high rotation rates. 

Work was begun on the development of improved second-order closure models. 
In all of the currently popular second-order closure models it is assumed that the 
deviatoric part of the velocity pressure-gradient correlation ollij and the deviatoric 
part of the dissipation rate corelation r»D ;; ar#> functions of r;j, Sij, and Clij. Hence, 
we considered the most general model of the form 

D^kl + D&kl = fkl(Tij^Sij,Qij) (11) 

Equation (8) should be subjected to the constraints: 

(i) Form invariance under a change of coordinates (c./., Smith 1971). 

(ii) Material frame-indifference in the limit of two-dimensional turbulence (see 
Speziale 1981). 

and the fact that D^iji and are traceless. This led to the most general form, 


£>IIfcl + D^kl —01~?{ T kl ~ + 02 kSlcl + 03 4 ( T km T ml 0 T mn T mn^kl) 

q J 3 q* 3 

~\~04 (Skm§ml ~ n^mn^mn^kl) 

( 3 2 

+/^5 Sjnl TimSmk ^^mn^mn^kl) 

6 2 

H o ( T km T mnSnl “I" T lm T mnS n k ~ ~T rnn T n pSp Tn fi kl ) 

<r 3 

+/?7 { T kmSmnSnl + T , | m 5 mn '5nfc — n T mnSnpSpm&kl) 

t O 

1 2 
+/?8 -~{Tkm T mnSnrS r l + T/ m T mn jS , n r*^rfc ‘Z^mn^np^prSrm^kl) 

€ 1 

+ 2(1 — f$$)(jkrrSllm + T lm^km) + ft 9 ~^( T fem T mn^nl + T lm T mn^nfc) 

q . (12) 

The coefficients 0i are functions of the invariants. This model is substantially 

simpler than previous attempts at general representations which contained several 
redundant terms (c.f., Reynolds 1987 for a summary of such previous representa- 
tions). In the limit of a two component turbulence, Realizability (Shih 1987, private 
communication) requires that 


01 = 2 - y^T kl T lk 0 3 

02 =0 

04 =0 

1 fc 1 

06 = (7 '^kl‘^lmSmk06 “I" TklSl m Smk07 f" — / Tkl^’lm^mnSnk06 )/fpqSpq 

k € € 


( 13 ) 
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and hence, in the first approximation this model has only five undetermined con- 
stants. This model will be investigated in the future in collaboration with Dr. T.-H. 
Shih. It is interesting to note that this model is consistent with the numerical re- 
sults of this study which predict that an initially isotropic turbulence in a rotating 
frame decays isotropically. 
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